Spectrum of the Dirac Hamiltonian with the mass-hedgehog in arbitrary dimension 
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It is shown that the square of the Dirac Hamiltonian with the isotropic mass-hedgehog potential 
in d dimensions is the number operator of fictitious bosons and fermions over d quantum states. 
This result allows one to obtain the complete spectrum and degeneracies of the Dirac Hamiltonian 
with the hedgehog mass configuration in any dimension. The result pertains to low-energy states 
in the core of a general superconducting or insulating vortex in graphene in two dimensions, and in 
the superconducting vortex at the topological - trivial insulator interface in three dimensions, for 
example. The spectrum in d = 2 is also understood in terms of the underlying accidental SU (2) 
symmetry and the supersymmetry of the Hamiltonian. 



I. INTRODUCTION 

The spectra of Dirac particle in topologically non- 
trivial backgrounds have been attracting attention ever 
since the pioneering work of Jackiw and Rebby on frac- 
tionalization of electric charge, [ij With the recent rise 
of graphene and topological insulators this problem has 
only gained in importance. The most interesting fea- 
ture of the spectrum probably is the appearance of 
states with precisely zero-energy, which are topologic ally 
protected and have several exotic manifestations. |2h13| 
The best understood example is provided by the two- 
dimensional Dirac equation in the vortex background, 
which in its various representations is applicable to both 
graphene ffllll - 13 and the surface of a topologi- 
cal insulator [3, [ij]- Very recently, it was shown [12.] that 
the unique zero-energy state also exists in the three- 
dimensional version of this Hamiltonian, in which the 
mass-vortex is replaced by its three-dimensional topo- 
logical equivalent, the hedgehog. This Hamiltonian de- 
scribes, for example, a superconductor coexisting with 
the interface between a topological and a trivial insula- 
tor, and its three-dimensional zero-modes were related to 
non-Abelian statistics for such topological defects. [T2|] 

Unlike the zero-modes, the rest of the energy spec- 
trum of topologically non-trivial Dirac Hamiltonians has 
seen relatively little study [l^, [3 • 1^ this paper we con- 
sider the spectrum of the Dirac Hamiltonian in the sim- 
plest topologically non-trivial background of a hedgehog 
in general number of d dimensions, in which the bound- 
ary of the configurational space wraps the Sd-i sphere of 
the order parameter (i. e. masses of Dirac fermion) space 
fully once. We show that the simplest Dirac Hamiltonian 
of this kind is in fact a square-root of the particle number 
operator for a collection of fictitious bosons and fermions 
distributed over d quantum states: 
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Here N = + the sum of the standard 

bosonic and fermionic number operators, p and r are the 
momentum and position operators, [fijPj] = iSij, and 
Ui and /3i, i — 1, ...d are the usual Hermitian Dirac ma- 
trices that satisfy the 2(i-dimensional Clifford algebra. 



{ai,aj} = {f3i,Pj} = 2Sij, {ai,/3j} = 0. This sur- 
prisingly simple formula, which may be understood as 
a generalization of the celebrated Dirac connection be- 
tween the bosonic number operator and the quantum 
harmonic oscillator, enables one to determine the spec- 
trum and the degeneracies of the Dirac Hamiltonian on 
the right hand side by purely combinatorial or algebraic 
means. In particular, it immediately follows that there 
is a non-degenerate zero-energy state in any dimension, 
which persists under smooth deformations of the mass- 
hedgehog configuration in Eq. (1). The above Hamil- 
tonian is particularly relevant in two dimensions, since 
there it represents the leading term in the expansion of 
any realistic Dirac vortex Hamiltonian near its center, 
and the spectrum that will be obtained here may be un- 
derstood as the first approximation to the one of the 
real system at energies close to zero. We therefore also 
provide an alternative understanding of the spectrum 
in terms of symmetries of the Hamiltonian, and discuss 
some consequences of these considerations. In particular, 
we show that the mathematical origin of the observed de- 
generacies are two separate accidental symmetries of the 
Hamiltonian: one closely related to the accidental sym- 
metry of the harmonic oscillator, and the second that 
may be understood as a form of supersymmetry. [l^ 



The paper is organized as follows. In the next section 
we first motivate and then define the hedgehog Dirac 
Hamiltonian. We then proceed to compute its spectrum 
in arbitrary dimension by introducing suitable bosonic 
and fermionic creation and annihilation operators. The 
anisotropic version of the Hamiltonian is solved in sec. 
HI. In sec. IV we discern a hidden SU{2) symmetry of 
the Hamiltonian, which explains a part of the observed 
degeneracies. In sec. V we present a set of operators 
which close the Clifford algebra, and which are respon- 
sible for the additional degeneracy of the spectrum. We 
close with the remark on applicability of our results to 
more realistic Hamiltonians that would represent vortices 
in systems with Dirac fermions, other comments, and the 
short summary. 
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II. BOSE AND FERMI OPERATORS AND THE 
SPECTRUM OF THE HEDGEHOG 



Inserting this into Eq. (6) the oscillator's zero-point en- 
ergy is canceled and 



To motivate what follows let us begin with the familiar 
Dirac Hamiltonian in a mass-vortex background in two 
dimensions [Mllil: 

Ha = c{aipi +a2P2) + |A(r)|(/3icos6'-f /32sin6'), (2) 

where and /3i, i — 1,2, are four-dimensional, unitary, 
Hermitian, anticommuting matrices, and (r, 0) are the 
usual polar coordinates. Finiteness of the energy of the 
vortex configuration, as usual, requires that the ampli- 
tude |A(r)| — ^ when r 0. In spirit of the harmonic 
approximation, expanding |A(r)| near the origin and re- 
taining only the leading term suggests an introduction of 
the Dirac hedgehog Hamiltonian as 



H = cd ■ p + mj3 ■ r, 



(3) 



in general number of spatial dimensions d. The last term 
represents a topologically non-trivial mass-hedgehog con- 
figuration in d dimensions, and in particular in d = 2, 
evidently H^ = H + 0{r'^). Here c is a velocity, and m 
is a parameter with the dimension of energy/length. We 
work in units in which 7i = 1, as usual, and, without loss 
of generality, will assume m > 0. Since in general the 2d 
Hermitian matrices ai and /3i form the Clifford algebra 
C{2d, 0), their smallest irreducible representation will be 
2'^-dimensional. 

Squaring the Hamiltonian H yields 



H^ = c^p + ni^f + imc/S ■ a. 



(4) 



Following Dirac, the bosonic operators may be defined as 
usual: 



k - (m/2c)i/2f, -Hi(c/2m)i/2p„ 



(5) 



so that — Sij, and the square of the Hamiltonian 

may be rewritten as 



H^ = TOc^(2S,|6, -t- 1 -I- i/3,a,). 



(6) 



To discern fermions in the problem one may recall that 
the creation and the annihilation operators for Majorana 
fermions [1] satisfy precisely the same Clifford algebra as 
the Dirac matrices. Therefore, if we define a set of new 
operators as 



di = (A -I- iai)/2, 



(7) 



the anticommutation rules for the Dirac matrices imply 
that these linear combinations act as fermionic annihila- 
tion and creation operators, since evidently {d^, dj} = 0, 
and {di, d]} = ^.y . Furthermore, the product of the Dirac 

matrices appearing in H^ can be simply written in terms 
of the same fermions as 



iPiai 



2d|di 



1. 



(8) 



H^ = 2toc^(6|&, -|-d|aO, 

i=l 



(9) 



as announced in the introduction for m ~ c ~ I. 

The spectrum of the hedgehog-Dirac Hamiltonian is 
therefore 



= ±V2Nmc 



(10) 

where N — OjjL 2... is an integer quantum number. This 
agrees with [ij where the spectrum was obtained by 
solving the differential equation in coordinate representa- 
tion. There is a non-degenerate state with precisely zero- 
energy, when the number of both bosons and fermions 
vanishes. Otherwise, the degeneracy of the eigenvalue 
E^ is 



Did,N) = lY^ 



dl 



{d + N -n - 1)! 



2^^id-n)M id~mN-ny. 



(11) 



The two factors under the sum represent the number of 
ways in which n fermions and N — n bosons can be dis- 
tributed over d quantum states, respectively. The factor 
of 1/2 in front originates in the "chiral" symmetry be- 
tween the positive and negative energy eigenstates, im- 
plied by the existence of the Hermitian operator 
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(12) 



which anticommutes with the Hamiltonian H. 

In the physical dimensions d — 1,2,3 for N ^ we 
this way find: 



D{l,N) = l, 
D{2,N) = 2N, 
D{3,N) = 2N^ + 1, 



(13) 
(14) 
(15) 



with the degeneracy in d dimensions being a polynomial 
function of N of the degree d — 1. In d=l all the levels 
are non-degenerate. For a general dimension one can also 
write, 



D{d,N) = 



{N + d-iy. 
1)! 



F{~d,~N,l~d~ N,-l) (16) 



2iV!(d- 

where F is the hypergeometric function 2F1. 

It also transpires that the eigenstates including the 
zero-mode can be chosen to be real. The Clifford algebra 
C{n,m) for n -\- m-odd has a real representation only if 
n — m = 1, [16] with the dimension 2''. The Hermitian 
matrices {ai, j3i,^}, J = 1, ■■■d can thus always be chosen 
so that d -|- 1 of them are purely real, and the remaining 
d purely imaginary. [l3| Let us therefore chose Ui to be 
imaginary and j3i to be real. The time-reversal operator 
that commutes with the real Hamiltonian in this rep- 
resentation is simply that of complex conjugation, and 
therefore all the eigenstates can be chosen to be real. 
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III. SPECTRUM OF THE ANISOTROPIC 
HEDGEHOG 

One can determine the spectrum of a more general 
Hamiltonian, used for example in ref. [l^ 

H = cd ■ ]} + P ■ n{r), (17) 

with Hi ~ MijVj, using a similar technique. We assumed 
here the velocity to be isotropic, but it is easy to see 
that even if that would not be the case a simple rescaling 
would map the problem onto the Hamiltonian in the last 
equation. Squaring the Hamiltonian now yields 

d 

= ^{c^P^P^ + n{M^M),jfj + icP.M.jaj). (18) 

i=l 

A g eneral matrix M can be written in the diadic form 
[III as 

d 

M = ^m,{vi® III), (19) 

i=l 

with — {fin, fi2i, ■■■fJ'di), i = 1, ...d as the orthonor- 
mal eigenbasis and {to^} the spectrum of the symmetric 
matrix M^M: 

d 

MTAf = ^m2(^,®^7). (20) 

i=l 

We may therefore still define the bosonic operators as 

h = {\m,\l2c)^/^h + i{c/2\m,\f'^p^, (21) 

and, for mi > 0, the fermionic as 

a, = 0, + ia,)/2, (22) 

where f; = J^kf^kifk, Pi = J^kJ^kiPk, di = Y^kt^kiO-k, 
and Pi = i^kiPk- If rrii < Q, d and /3 matrices need to 
be exchanged in the definition of the fermionic operator. 
In terms of these "particles" one finds that 

d 

= 2cY,\mi\{h% + a\a,). (23) 

Evidently, H still has a non-degenerate zero-mode, but 
the rest of the spectrum becomes modified into 

d 

E^ = ±{2cY,Nk\mk\)"\ (24) 

k=l 

where Nk is an integer, including zero. The degeneracy 
for a generic set of incommensurate {ruk} of these eigen- 
values is now at most 2"^^^. 



In terms of the particle operators the Hamiltonian it- 
self takes the form 

d 

H = (2c)i/2 ^ \m^'"\a% + h[ai), (25) 

i=l 

from which it is clear that it commutes with the total 
number of bosons and fermions in each state. To obtain 
the eigenvectors therefore one needs to diagonalize a ma- 
trix of typical size 2"^ x 2"^, or 2D{d,N) x 2D{d,N), in 
the isotropic case rrii = m. 

IV. ACCIDENTAL SU(2) SYMMETRY IN d = 2 

Let us now focus on the physically most obviously rel- 
evant two-dimensional case, d = 2, and understand the 
obtained degeneracies in the isotropic case in a different 
way. First, the complete vortex Hamiltonian Hq in Eq. 
(2) is invariant under rotations of the coordinate frame 
generated by 2J2, where 

J2 = L2 + ^2 (26) 

with L2 = r X p/2 and S2 — i{a2ai + /32/3i)/4. In terms 
of the particle operators, 

X2 = iiAci - cIc2)/2, (27) 

where c = b (boson) for X — L, and c — a (fermion) 
for X — S. Of course, the hedgehog Hamiltonian H 
is rotationally invariant as well. However, since H also 
commutes with the total particle number in each state 
it also commutes with their difference: [H, J3] — 0, with 
J3 = L3 + 5*3, and 

I3 = {c\ci - clc2)/2. (28) 

This is an accidental degeneracy, which originates in the 
linear dependence of the hedgehog Hamiltonian on coor- 
dinate and momentum. By Jacobi's identity H then also 
commutes with the commutator of the (exact) generator 
J2 and the (accidental) generator J3, which happens to 
be non-trivial: Ji = Li + Si, where 

Xi = {c\c2 + clci)/2. (29) 

Since, by our definition, 

2 

4<^5/3/2, (30) 

a, 13=1 

where ct' are Pauli matrices, it is easy to see that these 
operators close the algebra SU{2) x SU{2): 

[X„Yj] = ieijkXkSxY, (31) 
where X, y = L, S. 
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We can therefore diagonalize the Hamiltonian H , the 

Casirair operator J^, and one Jk together, as usuaL For 
convenience we will choose here the latter to be the acci- 
dental generator J3. (One can of course also choose the 
customary exact rotational generator J2, which becomes 
analogous to our choice after a certain redefinition of the 
particle operators [l^.) The source of degeneracy will 
now turn out to be twofold: first, there is the "usual" 
accidental degeneracy of the harmonic oscillator due to 
different allowed quantum numbers m (the eigenvalues 
of J3) for a given value of j (where the eigenvalues of 

J2 are j(j + !))• This is a consequence of the accidental 
SU{2) symmetry. Second, there is an extra degeneracy 
due to fermions, which can be understood as different al- 
lowed values of j at given energy. We will see in the next 
section that this degeneracy is implied by a hidden super- 
symmetry of the hedgehog Dirac Hamiltonian. For now, 
let us simply ask which values of the quantum number 
j can be present among the eigenstates with the energy 
= ^/2Nmc. It is useful again to consider rather 
than H, because of its larger symmetry: 

= = 0. (32) 

We can therefore label the eigenstates of with the 
quantum number TV with four additional quantum num- 
bers: {l,mi, s,ms}, in standard angular momentum no- 
tation. In this basis H itself is represented by a purely 
off-diagonal matrix. We first note that from the defini- 
tion of Li it readily follows that 

- (nj + n^)(nj + 7i\ + 2)/4, (33) 

with = 6|,6fe, and therefore we can recognize the quan- 
tum number I as being 

I = (34) 

where = N\ + N2 is the total number of bosons in 
a given state. Note that this reproduces the well-known 
result [III that the degeneracy of the E = {2N^ + 1) level 
oi d — 2 harmonic oscillator is 2Z -I- 1 = + 1. For the 
fermions, on the other hand, 

i2=3(^l_M)2. (35) 

We can therefore also write that s(s -I- 1) = 3mg, and as 
there are only two possible values of to^ in the problem, 
namely, and 1/4, in fact s — |ms|. Since this number 
is finite only when the total number of fermions — 
N[ + N2 = 1, and zero when = or = 2, we can 
finally write 

s = {2- Nf)/2. (36) 

Since for > 2, Z > s, and s = 0, 1/2 one can write 
all the allowed values of the the total quantum number j 
as 

j = l±s^{N-Nf ±Nf{2-Nf ))/2. (37) 



j^=l/2 ^ ' ' E/m j>=3/2 

— V?'- 

.)=l/2 
1 ■ 

,/=l/2 





FIG. 1: The bound state spectrum of hedgehog Dirac Hamil- 
tonian in d = 2. The energy is given by i5jv/cj = ±y/N where 
to = \/2mc. For the energy levels with > 1, there are two 
degenerate multiplets with the quantum number j — j> and 
j = j< given by Eqs. (38) and (39). 

In spite of the three possible values of the fermion number 
and two values of the sign in the last equation, for 
the total number of particles N > 2 there are in fact only 
two values of j allowed: 

J> = N/2, (38) 

j< = (iV/2) - 1. (39) 

The degeneracy of the energy level E^ — ±\^2N is there- 
fore 

Di2, N) = (2j> + 1) + (2j< + 1) = 2N. (40) 

For = 1, on the other hand, it is either that I — 0, 
s — 1/2, or I = 1/2, s — 0. In either case j = 1/2, so 
D{2, 1) = 2. This way the Eq. (14) becomes correct for 
all 7^ 0. The obtained spectrum and the degeneracies 
are represented in Fig. 1. 



V. HIDDEN SUPERSYMMETRY 

The appearance of two values of the quantum number 
j at the same energy is not explained by the accidental 
SU{2) symmetry from the last section. Although acci- 
dental, its role in the present problem is analogous to the 
role played by the exact rotational symmetry in the spec- 
trum of the hydrogen atom. The reader may recall that 
the appearance of different values of the angular quantum 
number in that case was a consequence of the existence 
of another conserved quantity in the problem, namely, 
of the Runge-Lenz-Pauli vector [ll]. We now show that 
a similar vector operator exists for the hedgehog Dirac 
Hamiltonian, and that it may be used to understand the 
degeneracies, and even derive the energy eigenvalues by 
purely algebraic means. 
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Let us now for simplicity set 2mc = 1, and define three 
additional Hermitian operators Aj, i = 1,2, 3, which also 
obey j4? = (without the summation convention over 
the repeated indices): 

2 

</3(«a^/3 + «/3^l)- (41) 

The following algebra of commutation and anticommu- 
tation relations is then readily found: 

{Ai,Aj} = 2HHij, (42) 
[Ai,Jj]=ieijkAk, (43) 

{H,Ai} = AJi. (44) 

It is relatively easy to sec that these relations then imply 
that [H,Ji] = 0, and [Jt,Jj] = icijkJk- The Eq. (42) 
says that within each eigenspace of H"^ the operators Ai 

close a Clifford algebra, and Eq. (43) says that A is in- 
deed a vector under the accidental SU{2) symmetry from 
the last section. Note also that the chirality operator 7 
satisfies, 

{j,Ai} = [j,Ji]=0. (45) 

We may then use the above algebra to demonstrate 
that in each energy subspace there could be at most two 
different values of the quantum number j, and that when 
both of the values coexist they will differ by unity. From 
Eq. (43) it follows that either 

A+\E^,j,m)^\E^,j',m+l), (46) 

or the left hand side is zero, where = Ai + iA^, 
and the first label of the state is the value of H"^, not of 
H. Therefore by choosing m = j we may immediately 
conclude that / > j + 1. Since from Eq. (43) we also 
have [J+, = 0, it follows that, if not zero, then 

A+\E^j,j}^\E^,j + l,j + l}, (47) 

i. e. in this case it must be that in fact )' = j + 1. 
Finally, since from Eq. (42) it is also true that = 0, 
it follows that in particular, 

Al\E',j,j)=0. (48) 

So in each eigenspace of there is either only one, or 
two allowed values of the quantum number j. Recalling 
that {7, H} = 0, from Eq. (45) it readily follows the 
same is true not only for , but for the Hamiltonian H 
itself. 

To derive the spectrum from the above algebra, con- 
sider the state j>, j>), where the first label is now the 



eigenvalue of the Hamiltonian H. Using Eq. (44) it is 
easy to see that 

{E,j>,j>\A3\E,j>,j>) = '^. (49) 

On the other hand since [A3, J3] — 0, the operator A3 
has finite matrix elements only between the states with 
the same eigenvalue of J3. If that eigenvalue is j> there 
are only two such states, and therefore in principle 

As\E,j>,j>) = ^\E,j>,j>) +x\- E,j>,j>), (50) 

where x is an unknown coefficient. On the other hand, 
Eq. (45) implies that 

{-E,j>,j>\A3\E,j>,j>) = -{E,j>,j>\A3\- E,j>,j>) 

(51) 

i. e. that the coefficient x is purely imaginary. Since on 

the other hand, the Hamiltonian H written in the form 
in Eq. (25) is manifestly real, the coefficient x must be 
real as well, and so we conclude that a; = 0. 

Finally, since A| = H^, we may write an equation for 
the allowed energy eigenvalues as 

= (^)^ (52) 

and recover the energy spectrum as E = ±y/2jy (in our 
current units 2mc =1), with the degeneracy 2)> + 1 + 
2(j> — 1) + 1 for j> > 1, and 2j> + 1 for j> < 1, where 
j> assumes all the allowed values, j> = 0, 1/2, 1, 3/2, 
as usual. 



VI. REMARKS 

If the mass order parameter is isotropic but it satu- 
rates, for example, ni(r) = mvi for r < L, and ni{r) = 
mL for r > L in Eq. (17), one can think of the spec- 
trum derived here as being approximately correct if the 
size of the defect's core L ^ y/c/m. This is because the 
low-energy eigenstates, with the energies ^ y^mc much 
below the continuum threshold at ^ mL, are Gaussian 
functions expon entially localized within the length scale 
of order ^/c/m, and as such will be little affected by the 
modification of the Hamiltonian far from the center of 
the vortex. We may thus expect to find D{d,N) levels 
for small N close together in energy. The zero-energy 
state, however, being non-degenerate, will still remain at 
zero in order to preserve the exact reflection symmetry 
between the positive and negative energy eigenstates of 
the Hamiltonian. 

Since in = 2 the Hamiltonian studied here can be 
understood as the leading term of a general Dirac vor- 
tex Hamiltonian near the origin, one may wonder how 
much of the degeneracy we identified will remain once 
the higher order terms in the expansion are included. 
Non-linear terms in coordinate will of course break the 
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accidental degeneracy, leaving only the U{1) rotational 
sub-symmetry. Since, as already discussed, by choosing 
imaginary a and real /3 the vortex Hamiltonian becomes 
real, the time-reversal operator is trivial, and by itself 
does not imply any residual degeneracy. We therefore 
expect the higher order terms to resolve the degeneracy 
completely. In the limit of large core where the linear 
approximation is adequate, the low-energy states should 
appear in bundles of 2N states grouped together. For 
a vortex in the superconducting (sl-flli ITsj . Neel Q or 
Kekule d, [2l|, order parameter in graphene, the pres- 
ence of spin degree of freedom of course doubles these 
near-degeneracies. In this case even the rotational sym- 
metry is only approximate (22j |. and the terms quadratic 
in momentum reduce it to C3. 

It may also be interesting to note the transformation 
properties of the generators of the SU{2) symmetry un- 
der time reversal. Since under time reversal the bosonic 
operators stay invariant, it is easiest to choose the repre- 
sentation of the Dirac matrices in which the same is true 
for the fermionic operators: this is again all /3i real and 
all ai imaginary. Time reversal in this representation is 
just the operation of complex conjugation, and whereas 
the generator J2 of real-space rotation is as usual odd, 
the accidental generators Ji and J3 are even. This trans- 
formation property is different from that of the true an- 
gular momentum, but still perfectly consistent with the 
enclosed SU{2) algebra. [Hj 



VII. CONCLUSION 

To summarize, we have found that the energy level of 
the linearized vortex Dirac Hamiltonian in two dimen- 
sions, = ±\/2N, is 2A^-degenerate for iV > 1. This 
degeneracy may be viewed as a consequence of an acci- 
dental SU(2) symmetry and a supersymetry, which to- 
gether lead to two multiplets, with the quantum numbers 
j = N/2 and j = {N/2) — 1, being degenerate in energy. 
We have also obtained the spectrum and the degenera- 
cies of the Hamiltonian in general dimension, showing in 
passing that there is always a non-degenerate zero mode. 

The bound-state spectrum obtained here may be ob- 
servable by scanning tunneling microscopy, with the su- 
perconductivity in graphene or on a surface of a topolog- 
ical insulator induced by a proximity effect, for example, 
as a finite local density of states well below the bulk gap. 
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